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^7^: ABSTRACT 
O 

^ ■ We studied the collapse of rotating molecular cloud cores with inclined mag- 
netic fields, based on three-dimensional numerical simulations. The numerical 



simulations start from a rotating Bonnor-Ebert isothermal cloud in a uniform 
magnetic field. The magnetic field is initially taken to be inclined from the ro- 
tation axis. As the cloud collapses, the magnetic field and rotation axis change 
their directions. When the rotation is slow and the magnetic field is relatively 
strong, the direction of the rotation axis changes to align with the magnetic field, 
as shown earlier by Matsumoto & Tomisaka. When the magnetic field is weak 
and the rotation is relatively fast, the magnetic field inclines to become perpen- 
dicular to the rotation axis. In other words, the evolution of the magnetic field 
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and rotation axis depends on the relative strength of the rotation and magnetic 
field. Magnetic braking acts to align the rotation axis and magnetic field, while 
the rotation causes the magnetic field to incline through dynamo action. The 
latter effect dominates the former when the ratio of the angular velocity to the 
magnetic field is larger than a critical value VLq/Bq > 0.39 G^/^ c^^, where Bq, 
flo, G, and denote the initial magnetic field, initial angular velocity, gravita- 
tional constant, and sound speed, respectively. When the rotation is relatively 
strong, the collapsing cloud forms a disk perpendicular to the rotation axis and 
the magnetic field becomes nearly parallel to the disk surface in the high den- 
sity region. A spiral structure appears due to the rotation and the wound-up 
magnetic field in the disk. 

Subject headings: ISM: clouds: ISM: magnetic fields — MHD — stars: formation 

1. Introduction 

Magnetic fields and rotation are believed to play important roles in the gravitational 
collapse of molecular cloud cores. For example, the outflow associated with a young star is 
believed to be related to the magnetic flcld and rotation of a protostar and its disk. The 
rotation of various molecular clouds have been studied by Casclli et al. (2002), Goodman 
et al. (1993), and Arquilla & Goldsmith (1986), who found that the rotation energy is not 
negligible compared with the gravitational energy. Crutcher et al. (1999) obtained strengths 
of the magnetic field for various molecular clouds by Zeeman splitting observations and also 
concluded that the magnetic energy of a molecular cloud is comparable to the gravitational 
energy. The direction of the magnetic field is also crucial for cloud evolution, because it 
controls the direction of outflow and the orientation of the disk. Polarization observations 
of young stellar objects suggest that circumstellar dust disks around young stars are aligned 
perpendicular to the magnetic fleld (e.g., Moneti 1984; Tamura & Sato 1989). However, 
the direction of a large-scale magnetic field of an ambient cloud and that of a small-scale 
magnetic field around a molecular core do not always coincide. For example, the Barnard 1 
cloud in Perseus exhibits field directions different from the ambient field in three of its four 
cores (Matthews & Wilson 2002). Recently, high-resolution observations by Matthews et al. 
(2005) have shown that the polarization angle measured in the OMC-3 region of the Orion 
A cloud changes systematically across the core. These findings indicate that the spatial 
configuration of magnetic field lines in a molecular core is not simple. 

Dorfi (1982, 1989) and Matsumoto & Tomisaka (2004, hereafter MT04) numerically 
investigated the contraction of a molecular cloud when the magnetic field lines are not 
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parallel to the rotation axis (non-aligned rotator). Dorfi (1982, 1989) showed that a formed 
disk changes its shape into a bar or a ring depending on the angle between the magnetic 
field and the rotation axis at the initial stage. However, the evolution was not calculated to 
high densities. MT04 reinvestigated the collapse of rotating molecular cloud cores threaded 
by oblique magnetic fields by high-resolution numerical simulations. They found that strong 
magnetic braking associated with outfiow causes the direction of the angular momentum to 
converge to that of the local magnetic field, resulting in convergence of the local magnetic 
field, angular momentum, outfiow, and disk orientation. However, their study was restricted 
to clouds with an intermediate rotation rate {Q — 7.11 x 10~^ yr~^). It is expected that the 
evolution is affected not only by the magnetic field strength but also by the rotation rate. 
In this paper, we extend the parameter range of the MT04 study and explore the evolution 
of rotating magnetized clouds more generally. 

In the case of a cloud in which the magnetic field is parallel to the rotation axis (aligned 
rotator), four distinct evolutions are observed according to the magnetic field strength and 
rotation rate (Machida ct al. 2004, 2005a,b, hereafter Papers I, II, and HI). In the isothermal 
regime a contracting disk is formed perpendicular to the magnetic field and the rotation axis. 
In the disk, the magnetic field strength, angular rotation speed, and gas density are correlated 
with one another and satisfy the magnetic flux-spin relation: 

(0.2)2 x'47rGp, + (0.36)2 x'gTrc^p, = ' ^' 

where Vtzci B^c, pci Cs, and G arc the angular velocity, magnetic flux density, gas density at 
the center, isothermal sound speed, and gravitational constant, respectively. 

In the case of a weak magnetic field and a small rotation rate, -F < 1, the cloud contracts 
spherically (spherical collapse). This increases F and a self-similarly contracting disk forms 
when F ~ 1 is reached. On the other hand, a cloud with a strong magnetic field and/or a 
fast rotation rate, F > 1, contracts only in the direction of the magnetic field and rotation 
axis (vertical collapse). This reduces F and a self-similarly contracting disk forms when 
F ~ 1. Hence, F controls the mode of contraction (spherical or vertical collapse). This is 
understood as follows: for F < 1 the support forces are deficient (support-deficient model), 
while for F > 1 the cloud is supported laterally by rotation and/or the magnetic field 
(support-sufficient model) . 

The evolution can be further divided into two categories depending on whether the 
magnetic or centrifugal forces are more effective in forming the disk. A model with a spin- 
to- magnetic fiux ratio Qq/Bq larger than a critical value {Q,/B)crit — 0.39G^^'^cJ^ forms a 
disk mainly supported by the centrifugal force, while that with a smaller ratio forms a disk 
by the Lorentz force. 
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In addition to the differences in the evolution of the pre-protostellar phase (isothermal 
phase), subsequent evolution is affected by the magnetic field and rotation. A first core 
consisting of adiabatic H2 molecular gas experiences fragmentation only if the initial cloud 
is rotation dominated, (Qq/Bq) > {Q/B)c.rit- In the support-deficient regime (F < 1), 
fragmentation proceeds through a deformation forming a ring, while in the support-sufficient 
regime {F > 1) fragmentation from a bar appears as well as ring fragmentation. 

Based upon our previous results (Papers 11 and 111), all models of Dorfi (1982, 1989) 
and MT04 belong to a type of "support-sufficient", F > 1, and "magnetic-force-dominant", 
Qq/Bo < {Q/B)ait, models. In this paper, wc investigate the evolution of a rotating isother- 
mal cloud with an inclined magnetic field, not restricted to the above type, in large parameter 
space, and also explore three other models which have not been previously studied. 

The plan of the paper is as follows. The numerical method of our computations and the 
framework of our models are given in §2 and the numerical results are presented in §3. We 
discuss the geometry of the collapse in §4, and compare our results with previous works and 
observations in §5. 



2. Numerical Model 



Our initial settings are almost the same as those of MT04. To study the cloud evolution 
and disk formation, we use the three-dimensional magnetohydrodynamical (MHD) nested 
grid method. We assume ideal MHD equations including the self-gravity: 



| + V.(^) = 0. 

d 1 

+ p{v ■ V)v = -VP - —B X (V X S) - pV0, 



dt 



dB 

'dt 



An 

V X (v X B), 



(2) 

(3) 

(4) 
(5) 



where p, v, P, B, and denote the density, velocity, pressure, magnetic flux density, and 
gravitational potential, respectively. The gas pressure is assumed to be expressed by the gas 
density (barotropic gas) as 

^ Peri , 



1 + 



2/5' 



(6) 



where = 190ms ^ and pc 



1.9205 X 10-1=^ g cm-=^ {n^ 



5 X 10^° cm ^ for an assumed 



mean molecular weight of 2.3). This equation of state implies that the gas is isothermal at 
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T = 10 K for n <C ricri and is adiabatic for n » ncri (Masunaga & Inutsuka 2000). For 
convenience, we define the core formation epoch as that for which the central density {ric) 
exceeds ricri- We also call the period for which ric < ricri the isothermal phase, and the period 
for which ric ^ ^cri the adiabatic phase. 

In this paper, a spherical cloud with a critical Bonnor-Ebert (Ebert 1955; Bonnor 1956) 
density profile, pbei is assumed as the initial condition, although a filamentary cloud is 
assumed in Papers l-lll. The cloud rotates rigidly (Qq) around the ^-axis and has a imiform 
magnetic field (Bq). To promote contraction, we increase the density by a factor / (density 
enhancement factor) as 

. . ^ r PBE(r) / for r < R^, 
I PBE(Rc)f for r>Rc, 

where r and Rc denote the radius and the critical radius for a Bonnor-Ebert sphere, re- 
spectively. We assume Pbe(O) — 7.39 x 10~^^g cm~^, which corresponds to a central num- 
ber density of nc,o — 5 x 10^ cm"^. Thus, the critical radius for a Bonnor-Ebert sphere 
i?c = 6.45c,[47rGpBE(0)]-V2 corresponds to Rc = 2.05 x lO'' AU for our settings. 

The initial model is characterized by three nondimensional parameters: a, cu, and ^o- 
The magnetic field strength and rotation rate are scaled using a central density po — Pbe(O)/ 
as 

a^B^J{Anpocl), (8) 

and 

a; = Qo/(47rG'po)'/'- (9) 

The parameter 6*0 represents the angle between the magnetic field and the rotation axis 
(2;-axis). Thus, the initial magnetic field is given by 



(10) 



in Cartesian coordinates {x, y, z). The above definitions of a and cu are the same as those 
of Papers I-III. Although we added a finite non-axisymmetric perturbation in Papers I-III, 
we assume no explicit non-axisymmetric perturbation in this paper. 

We calculated 36 models, widely covering the parameter space. 20 typical models are 
listed in Table 1. The model parameters (a, a;, and ^o); density enhancement factor (/); 
ratio of the thermal {ao), rotational (/?o), and magnetic (70) energies to the gravitational 
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energy;^ initial central number density (no), magnetic field strength (Bq); angular velocity 
(Qq); and total mass inside the critical radius (r < R^) arc summarized in this table. The 
models SF, MF, and WF in MT04 are also listed. The clouds in MT04 have stronger (models 
SF and MF) or equivalent (model WF) magnetic fields compared to our models and have 
an intermediate rotation rate. 

We calculated cloud evolutions up to n ~ 10^"^cm~^ using an ideal MHD approxima- 
tion. The ideal MHD approximation is fairly good as long as the gas density is lower than 
~ 10^^ cm^^ (Nakano 1976; Nakano et al. 2002). However, ohmic dissipation affects proto- 
stcllar collapse, especially at high densities exceeding n ~ lO^^cm"^ (Nakano et al. 2002). 
Nakano et al. (2002) have shown that the magnetic field is coupled with gas for n < 10^^ cm~^, 
indicating that the assumption of an ideal MHD is valid in the isothermal phase. In the adia- 
batic phase, the number density in the adiabatic core exceeds ~ 10^^ cm~^, and the magnetic 
field begins to decouple as the ohmic dissipation being effective. Our simulation may there- 
fore overestimate the magnetic field, especially for a dense core (n ^ lO^^cm"^). Since 
we are interested in the direction of the magnetic field, rotation axis, and disk normal in 
the isothermal phase, we show the results of cloud evolution mainly for low density cores 
{n <^ 10^^ cm~^) in which the ideal MHD approximation is valid. The effect of ohmic dissi- 
pation and cloud evolution, outfiow, and jets in high density cores will be investigated in a 
subsequent paper. 

We adopt the nested grid method (for details, see Paper 11). The nested grid consists 
of concentric hierarchical rectangular grids to give a high spatial resolution near the origin. 
Each level of the rectangular grid has the same number of cells (128 x 128 x 128), but the 
cell width h{l) depends on the grid level I. The cell width is reduced by a factor of 1/2 as the 
grid level increases by 1 {I ^ I + 1). We begin our calculations with 6 grid levels {I — 1-6). 
The box size of the initial finest grid I — 6 was chosen to be 2Rc. The coarsest grid {I — 1), 
therefore, has a box size equal to 2^ R^. A boundary condition is imposed at r = 2^ Rf., 
where the magnetic field and ambient gas rotate at an angular velocity of Qq (for details see 
MT04). Due to this large simulation box, it takes t ~ 40 free-fall time until the Alfven wave 
generated at the cloud center reaches the simulation boundary, even in the model with the 
strongest magnetic field. Hence, the boundary condition does not affect the central cloud 
because our calculations end within ~ 10 free-fall time. The highest level of grids changes 
dynamically. A new finer grid was generated whenever the minimum local Jeans length Aj 
become smaller than 8/i(Zmax), where h is the cell width. The maximum level of grids was 
restricted to Zmax = 20 in typical models. Since the density is highest in the finest grid, the 



^Denoting the thermal, rotational, magnetic, and gravitational energies as U, K, M, and W, the relative 
factors against the gravitational energy are defined as ao = U/\W\, /3o = -fC/IW^I, and 70 = M/|W^|. 
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generation of new grid ensures the Jeans condition of Truelove et al. (1997) with a margin 
of a factor of 2. We adopted the hyperbohc divergence B cleaning method of Dedner et al. 
(2002). 

3. Results 

Our models are characterized mainly by the strength of the magnetic field (a) and the 
angular velocity (u). We calculated six groups of models, groups A-F, distinguished by the 

values of a and u (see Table 1). The models are designated by group (A, B, C, D, E, or F) 
and ^0 (0°, 30°, 45°, or 60°). Hence, model AOO has a = 0.01, u = 0.01, and 60 = 0°, while 
model A45 has the same a and u, but ^0 = 45°. Models with ^0 = 0° ^''^^ "aligned rotators" 
and those with 6q 7^ 0° are "non-aligned rotators." The aligned rotator models of groups 
A, B, C, and D (AOO, BOO, COO, and DOO) have the same a and uj as groups A, B, C, and 
D in Papers II and III and have the same magnetic field strength and angular velocity at 
the center, although the distributions of the density, magnetic field, and angular velocity are 
different. Groups E and F are newly added in this paper. 

The evolution of the aligned rotator models can be divided into four types (§1) by two 
criteria: F ^ 1 and VLq/Bq ^ 0.39G^/^cJ^. Groups A, B, and F are magnetic-force-dominant, 
Q0/-60 < 0.396*^/^07^ ^ = 0, and Groups C, D, and E are rotation-dominant. The models 
are described in the following subsections (§3.1-§3.3). 

3.1. Magnetic-Force-Dominant Disks 

In this subsection, we show the evolutions of groups A, B, and F. These groups are 
expected to form magnetic-force-dominant disks in the isothermal phase since Q^q/Bq < 
0.39G^/^cJ^. Before disk formation, the cloud is expected to collapse spherically in groups A 
and B, since they are support-deficient models, F ^ 1. On the other hand, one-dimensional 
collapse along the magnetic field lines is expected for model F, since this is a support-sufficient 
model, -F ;^ 1. 

3.1.1. Support- deficient models 

In this section, we consider the evolution of groups A and B. Figures 1 (a)-(d) shows 
the evolution for model A45. Model A45 has the parameters a — 0.01, u; — 0.01, and 
9o — 45°. This cloud rotates slowly {u; — 0.01) around the 2;-axis, and has a weak magnetic 
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field {a — 0.01). Figure 1 (a) shows the initial state of model A45. The cloud is threaded by 
a magnetic field running in the direction 9q = 45°. Figure 1 (b) shows the structure when 
the central density reaches Uc = 5.9 x 10^ cm"'^. Figure 1 (b) (/ = 9) covers (1/8)^ of the 
volume of Figure 1 (a) (/ = 6). Figure 1 (b) shows that the magnetic field lines run in a 
direction at an angle 6b — 45° from the z-axis, similar to the large-scale view (/ = 6) of 
Figure 1 (a), though they are squeezed near the center. The green colored disks in Figures 1 
(b), (c), and (d) indicate regions of p > (1/100) Pc on the mid-plane parallel to the disk-like 
structure (perpendicular to the disk normal). From the density contour lines in the x-z 
plane of Figure 1 (b) , it can be seen that the high density region is shghtly flattened and the 
cloud collapses along the magnetic held lines. We derive three principal axes (ai > 02 > as) 
from the moment of inertia according to MT04. We define the shortest axis as to be the 
disk normal axis and define ai and 02 to be the long and short axes on the disk mid-plane. 
The oblateness and axis-ratio of the high-density region of p > O.lpc are defined here as 
£oh = (0102) ^''^/os and £ar = 01/02 — 1, respectively. 

Figure 2 plots the oblateness (upper panel) and axis-ratio (lower panel) against the 
central density for group A. For model A45, the oblateness reaches Sob — 1-1 at the epoch 
of Figure 1 (b) {ric = 5.9 x lO^cm"^). This means that the high-density region maintains a 
spherical structure at this epoch. The axis-ratio grows only to ^ar — 10^'^ at the same epoch. 

Figure 1 (c) shows the high-density region at the core formation epoch (n^ = 5 x 
10^° cm^'^). The angle between the disk normal and the 2;-axis has increased to 6p ~ 41°, 
although that between the magnetic field lines and the 2;-axis is 9b — 45°, similar to Figure 1 
(b). Figure 3 shows the loci of the magnetic field B {9b,4>b), the rotation axis f2 {9^,(f>n), 
and the normal vector of the disk P {9p, (pp), all averaged within p > O.lpc- Each vector is 
projected on the z — plane and the distance from the origin represents the angle between 
a given vector and the 2;-axis. For example, a vector parallel to the y-axis is plotted at 
(0°, 90°). The dotted fine in Figure 3 indicates that the magnetic field rotates around the 
2;- axis from 0^ = 0° to 125° during the calculation, keeping the same angle with respect 
to the z-axis {9b — 45°). It can also be seen that the disk normal (the solid line) stays 
parallel to the magnetic field after ric > 10^ cm^'^, showing that a disk structure is formed 
perpendicular to the magnetic field and rotates with the magnetic field. On the other hand, 
the rotation axis, represented by the broken line, points in the direction 9q ^ 15° before 
ric ^ 5 X 10^° cm~^. Thus, the rotation axis maintains its initial direction in the isothermal 
phase. We summarize the angles of the magnetic field {9b, (Pb), rotation axis {9a, 0n), and 
disk normal {9p, (j)p) at the end of the isothermal phase in Table 2. The angles between the 
magnetic field and the rotation axis V'Bn, those between the magnetic field and disk normal 
■0BP, and those between the rotation axis and the disk normal -^qp are also listed in Table 2. 
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For example, the angle ■0bq is calculated as 

_ABxn\ 

^Bn = sm -J^^^- (11) 

We can confirm from Table 2 that the magnetic field and disk normal have almost the same 
directions {ipBP = 4°) and that the disk normal is inclined from the rotation axis at an angle 
of ipnp = 42° at the end of the isothermal phase. The oblateness reaches ^ob = 1-45 at the 
core formation epoch (Figure 2 upper panel). This oblateness is smaller than that of model 
AS (sob = 2.9) of Paper II. The increase in the initial central density from 5 x 10^ cm~^ (Paper 
II) to 5 X 10^ cm~^ suppresses the growth of oblateness. The disk forms slowly compared with 
models with large a or uj, similar to Paper II. The axis-ratio grows only to £ar = 7.1 x 10~^ 
at the end of the isothermal phase. 

Figure 1 (d) shows the central region at 140 yr after the core formation epoch (ric — ricn)- 
At this stage, the directions of the magnetic field and disk normal are parallel to each other 
and co-rotate around the 2;-axis, keeping the angles 9b, dp — 45°. The magnetic field and disk 
normal vector continue to rotate up to (f)B,(f)p — 125° at the end of the calculation period 
(n^ ~ lO"*^^ cm^^). The rotation axis begins to move away from the 2;- axis at the beginning of 
the adiabatic phase, and has an angle 6q ~ 20° with respect to the 2;-axis at the end of the 
calculation (asterisk). It can be seen that the three axes {B, O, and P) are in alignment in 
the adiabatic stage, as seen in MT04. Although the disk becomes thinner in the adiabatic 
phase, the oblateness is no more than £ob — 2, as shown in Figure 2. When the central density 
exceeds ric ^ 10^^ cm~^, the oblateness gradually decreases and the central region becomes 
spherical owing to the thermal pressure. Figure 2 indicates that the oblateness depends only 
slightly on Oq. On the other hand, the evolution of the axis-ratio clearly depends on ^o- 
The axis-ratio grows faster in models with larger 6q. The growth of the axis-ratio must be 
due to non-axisymmetry in the plane perpendicular to the magnetic field in magnetic-force- 
dominated models. Such non-axisymmetry is induced by the centrifugal force due to the 
rotation motion with O in the disk plane, which vanishes in an aligned rotator for which 
9o — 0. Thus, the axis-ratio does not grow in model AOO, for which the initial magnetic 
field lines are parallel to the rotation axis. This tendency is more marked in group C. We 
will discuss the correlation of the axis-ratio and the initial angle ^0 in §3.2. Although the 
high-density region rotates slightly around the rotation axis Q in the isothermal collapse 
phase [Figure 1 (a)-(c)], it begins to rotate in the adiabatic phase [Figure 1 (c) and (d)] 
when the gravitational collapse is suppressed by the thermal pressure. Comparing Figures 1 
(c) and (d), we can see the magnetic field lines and the disk-normal rotate around the z-a.xis. 

Figure 4 shows the evolution of the angles 9b, 9q,, and 9p against the central density for 
models AOO, A30, A45, and A60. The angles 0b, and 0p for model A45 are also plotted. 
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in this figure, we cannot see any differences in tfie angles Ob, Oq, and 9p of model AOO, because 
these angles are all zero. The angles 6b (broken lines) show that the magnetic field hardly 
changes with respect to the z-axis until the end of the calculation period. Thus, the magnetic 
field maintains its initial zenithal angle for group A. In contrast, the disk normal changes to 
the direction of the magnetic field after the central density exceeds n ^ 10^ — 10^ cm^"^ for 
all models AOO, A30, A45, and A60. Hence, we can see that the disk is formed perpendicular 
to the magnetic field, irrespective of the initial angle Oq. On the other hand, the rotation 
axis maintains its initial direction (z-axis) in the isothermal phase and then begins to move 
away from the z-axis in the adiabatic phase. The angles 0b and (pp increase slightly in the 
isothermal phase because the clouds rotate slowly. Figures 2 and 4 clearly show that the 
evolutions of the oblateness and the angles 9b, 9q, and do not qualitatively depend on 
the initial angle while the angles of the disk normal do depend on 9q. 

At the end of the calculation period for model A45, a disk structure is found in the 

region r ^ 200 AU, where the cloud has an oblateness of ^ob > 1-5. In this region, the 
direction of the disk normal varies depending on the disk scale. The disk normal is directed 
toward 9p ~ 45°, which is almost parallel to the magnetic field for r ^ 50 AU. On the other 
hand, for a disk in the range r ~ 100 - 200 AU, 9p increases from ~ 55° (r <^ 100 AU) 
to ~ 60° (r ^ 200 AU). Thus, the disk normal gradually becomes inclined and approaches 
the magnetic field direction, moving towards the center. The magnetic field strength also 
depends on the scale. It increases as the center is approached and has a maximum at the 
center. 

Next, we focus on group B (a = 0.1, = 0.01) which has a magnetic field 10^/^-times 
stronger than group A. The left panel of Figure 5 shows the same information as Figure 3 
for model B45 (^o — 45°) . The cloud structure, magnetic field lines, and velocity vectors for 
model B45 at the end of the isothermal phase are shown in the right panel, for which the 
contour, streamline, isosurface, and other notations have the same meanings as in Figure 1. 
In group B, a disk is formed perpendicular to the magnetic field, similar to group A. The 
direction of the magnetic field rotates around the ^-axis, keeping the initial angle 9^, similar to 
model A45. The directions of the magnetic field and the disk normal also coincide for model 
B45. The angles (pB and 0p, however, increase slightly in the isothermal phase, compared 
with model A45. Since the initial rotation speed uj — 0.01 is common for groups A and B, 
this difference must be due to the growth rate of the angular velocity fl, which is smaller in 
group B than in group A, as shown in Paper II. The growth rates of the angular velocity (Q) 
and magnetic fiux density (B) are large when the cloud collapses spherically (oc p^^^), while 
they are small when the cloud collapses laterally (oc p^^^) (Paper II). A cloud with a strong 
magnetic field (group B) forms a disk earlier than that with a weak magnetic field (group 
A), and thus, the growth rate of the angular velocity in group B decreases from fl oc p^/^ to 
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Q oc p^/^ at an earlier epoch than for group A. Further, magnetic braking is more effective in 
group B than in group A, since the initial magnetic field is stronger in group B. Model B45 
has 6q = 7° at the end of isothermal phase, while model A45 has 6^ = 2° (Table 2). That 
is, the rotation axis is more inclined with respect to the 2;-axis in model B45 than in model 
A45. This inclination is caused by magnetic braking (MT04). MT04 show that magnetic 
braking works more effectively for the perpendicular component of angular momentum to 
the magnetic field than for the parallel component. The angular velocity in model B45 is 
slower than model A45 for the following two reasons, prompt disk formation and effective 
magnetic braking. The direction of the rotation axis oscillates violently around the z-axis in 
the adiabatic phase (Figure 5) because magnetic braking is so effective, as noted in MT04 
(model MF) and Paper II. In conclusion, groups A and B exhibit a disk perpendicular to 
the magnetic field line irrespective of the initial angle ^o- 



3.1.2. Support- sufficient models 

Group F models are magnetic- force dominant, similar to groups A and B but with a 
strong magnetic field. Group F has the parameters a — 1 and uj — 0.05. As group F 
clouds have a stronger magnetic field and rapid rotation at the initial stage, to promote 
cloud contraction we set a larger density enhancement factor for group F (/ = 5) compared 
with that of groups A and B (/ = 1.68). We can confirm that the density enhancement 
factor does not greatly affect cloud evolution (see models A45 and [n], or models B45 and 
[o] in Table 2). The models in group F form disks through the magnetic force, as for models 
A45 and B45. For group F models, the cloud collapses along the magnetic field fines and 
contraction crossing the magnetic field lines is suppressed by the strong Lorentz force, while 
for models A45 and B45, the cloud collapses spherically. The locus of the disk normal (P) 
traces that of the magnetic field [B) in the left panel of Figure 6. The rotation axis is once 
inclined to 9^, ~ 70° in the isothermal phase, and then reverts to Oq, = 10° at the end of the 
isothermal phase. The loci of P, and ft are similar to those of models MF70 and MF80 
of MT04, which have parameters {a, u) = (0.76, 0.14). Although the three vectors {B, P, 
and O) do not completely converge to align in models MF70 and MF80, these vectors are 
roughly parallel, as seen in the left panel of Figure 6. The right panel of Figure 6 shows the 
cloud structure, magnetic field lines, and velocity vectors at the end of the isothermal phase. 
This panel shows that a disk forms perpendicular to the magnetic field lines, as for groups 
A and B. 

The axis-ratio grows to £ar = 0.4 in model F45 at the end of the isothermal phase 
(Table 2). The non-axisymmetric structure is caused by the centrifugal force whose rotation 
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vector is perpendicular to the magnetic field in magnetic-force-dominant models. The axis- 
ratio begins to grow after a thin disk is formed. The axis-ratios in group F are larger than 
those of groups A and B (see Table 2) as the disk formation epoch in group F is earlier than 
that of groups A and B. 

In groups A, B, and F (magnetic-force-dominant models), a disk is formed perpendicular 
to the local magnetic field. Our results agree quahtatively with those of MT04 (models SF, 
MF, and WF). 

3.2. Rotation-Dominant Disks 

In this subsection, we show the cloud evolution of groups C and E. These groups are 
expected to form rotation-dominant disks in the isothermal phase since Qq/Bq > O.SdG^^'^cJ^. 
The cloud collapses along the rotation axis (vertical collapse) in group C (support-sufficient 
models; F ;^ 1), while it collapses spherically in group E (support-deficient models; F 1). 

3.2.1. Support- deficient models 

The models of group E are rotation dominant (Qq/Bq > O.SQG^/^cJ^), although they 
have a slow rotation rate of a; = 0.05. Compared with group A, group E has a magnetic 
field 10^/^-times smaller but a 5-times larger angular velocity. Model E45 has parameters 
a = 0.001, uj = 0.05, and 9o = 45°. 

The models of group E form a disk through the centrifugal force. The left panel of 
Figure 7 shows that the rotation axis maintains its initial direction 6fi ~ 0° and the disk 
normal is also parallel to the 2;-axis (the rotation axis) in the isothermal phase. Therefore 
a disk forms perpendicular to the rotation axis, not along the magnetic field. The angle 6b 
increases gradually in the isothermal phase. This means that the magnetic field tends to be 
aligned along the direction perpendicular to the rotation axis. However, the isothermal phase 
ends before the magnetic field lines are completely directed in the perpendicular direction. 
Since the models in group E have initial states inside the B-Q relation hne,^ the magnetic 
field gradually becomes inclined as the cloud collapses slowly in a spherically symmetric 
fashion. Although the direction of the magnetic field does not coincide with the direction of 



^The cloud collapses slowly in a spherically symmetric fashion inside the B-fl relation line. Otherwise, 
the cloud rapidly collapses along the vertical axis when the model is distributed outside the B-Q relation 
line. 
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rotation and of the disk normal in the isothermal phase, the magnetic field, rotation axis, and 
disk normal begin to converge with each other after the core formation epoch. As shown in 
MT04, magnetic braking works preferentially for the component of the angular momentum 
perpendicular to the magnetic field, which drives the alignment of B and fl. 

3.2.2. Support- sufficient models 

Figures 8 (a)-(d) show the cloud evolution of model C45 in views along the y-axis (edge- 
on view; upper panels) and along the 2;-axis (face-on view; lower panels). Model C45 has 
the parameters a = 0.01, u = 0.5, and Oq = 45°. Compared with group A, group C has the 
same magnetic field but has a 50-times larger angular velocity. 

Figure 8 (a) shows the cloud structure at ric — 6.3 x lO^ cm"^. The spherical cloud 
collapses along the rotation axis (i.e. the z-axis), and then an oblate core forms at the center 
(upper panel). The magnetic field lines are slightly squeezed at the center (upper panel) and 
are rotated (f)B — 45° from the initial stage 0b = 0° [Figure 8 (a) lower panel]. Figure 8 (b) 
shows the core shape at ric = Q x 10^ cm~^. A thin disk is formed in the x-y plane. In this 
model, an extremely thin disk forms in the isothermal phase (^ob — 10 at rZc = 9 x 10^ cm~^, 
seen in Figure 9). In model CS of Paper II, a thin disk forms promptly in the early isothermal 
phase, because the lateral collapse is suppressed by a strong centrifugal force and hence the 
cloud collapses only vertically along the 2;-axis. In model C45, the cloud similarly collapses 
vertically and forms a disk promptly in the isothermal phase. The weak magnetic field in 
this model hardly affects the cloud evolution. Moreover, the magnetic field is compressed 
in the direction of the rotation axis and begins to run along the disk surface, as shown in 
Figure 8 (b). Outside the thin disk, the magnetic field lines run in the direction 6b — 45°. 
That is, the magnetic field lines emerge from the lower-left side of the cloud and escape from 
the disk in the upper-right direction. This configuration of the magnetic field appears only 
in non-aligned rotators. In contrast, the disk is vertically threaded by the magnetic field 
along the rotation axis in model COO. This configuration of the magnetic field is seen in all 
the models studied in Paper II, which was restricted to aligned rotators. 

Figure 10 shows the direction oi B, P, and Q for model C45. The inset at the lower-left 
corner is an enlarged view of the center. This shows that the direction of the magnetic field 
gradually moves away from the z-axis and towards 9b — 90°. Then, the direction of the 
magnetic field rotates around the ^-axis, keeping an angle of — 90°. On the other hand, 
the rotation axis hardly changes its direction from the initial state and remains directed 
along the z-axis. The disk normal is also oriented along the 2;-axis (i.e. the rotation axis). 
From Figure 8 (b) , it can be seen that a disk forms by the effect of the rotation and the disk 
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normal direction coincides with the rotation axis. 

Figure 8 (c) shows the central region at the core formation epoch (ric — 2.3 x 10^^ cm~^). 
It can be seen from this figure that a non-axisymmetric structure has formed and the central 
core has changed its shape from a circular disk [lower panels of Figures 8 (a) and (b)] to a bar 
[Figure 8 (c) lower panel]. The magnetic field lines run laterally, i.e. \Br\, jS^I ^ \Bz\, in the 
adiabatic phase [Figure 8 (c) and (d)]. Figure 8 (d) shows an adiabatic core when the central 
density has reached Uc = 6.9 x 10^^ cm^"^. A spiral structure is seen in this figure, which 
indicates that a non-axisymmetric pattern has formed, even if no explicit non-axisymmetric 
perturbation is assumed at the initial stage. (Although the non-axisymmetric patterns also 
appear in some models of Papers I-III, it should be noted that these patterns are due to a 
non-axisymmetric perturbation added to the density and magnetic field at the initial stage.) 
The magnetic field lines are considerably twisted in Figure 8 (d) . It should be noted that in 
this model, the inclined magnetic field induces non-axisymmetric perturbations, on behalf 
of the initial explicit perturbation. 

Figure 11 shows the magnetic field lines, the shape of the core, and the velocity vectors 
on the z — plane in the adiabatic phase for model COO. This figure shows that a ring is 
formed, as found in Paper III, without any growth of a non-axisymmetric pattern. In Papers 
I-III, we assumed a cylindrical cloud in hydrostatic equilibrium, in which the magnetic field 
and angular velocity arc functions of the radius r in cylindrical coordinates. On the other 
hand, the cloud is assumed to be spherical with a uniform magnetic field and angular velocity 
at the initial stage in model COO. In spite of these differences, a similar ring structure appears 
in both models COO and CS of Paper II. The lower panel of Figure 9 plots the evolution of 
the axis-ratio against the central density for group C. The axis-ratios for models C30, C45, 
and C60 begin to grow after a thin disk is formed (ric 5 x 10^ cm~^) and reach Egj ~ 0.5 
at the core formation epoch. The axis-ratio grows to £ar — 1 at ric — 10^^ cm~^ in models 
C30, C45, and C60, while no non-axisymmetric pattern appears in model COO. This shows 
that the non-axisymmetric pattern arises from the anisotropy of the Lorentz force around the 
rotation axis. A bar structure is formed by the non-axisymmetric force exerted by the inclined 
magnetic field, as shown in Figures 8 (b)-(d). This is confirmed by the fact that the short 
axis of the density distribution on the z = plane (the disk mid-plane) and the bar pattern 
rotate together with the magnetic field fines. The axis-ratio (the non-axisymmetricity) grows 
in proportion to p^^^ (10^ ^ ric ^ 10^°cm~^ in the lower panel of Figure 9), as found by 
Hanawa & Matsumoto (1999). Since the lateral component of the magnetic field {\B\ sin^o) 
is large (Figure 9 lower panel), the axis-ratio grows more in models with large Oq. 

The evolution of the angles 9b, On, dp, and (pB for group C are plotted against the central 
density in Figure 12. The angle between the magnetic field and 2;-axis becomes 9b — 90° 
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even in the early phase of isothermal collapse for all the models C30, C45, and C60. The 
rotation axis and the disk normal maintain their angles 6^, 9p ~ 0°. Figures 4 and 12 show 
that in both magnetic- and rotation-dominant models the directions of the magnetic field, 
rotation axis, and disk normal are qualitatively the same for models with the same a and a;, 
irrespective of 6*0 in the range 30° <9q< 60°. 

Figure 13 shows the magnetic field lines, velocity vectors, and density distribution for 
the epoch t = 1.52 x 10^ yr {ric = 1.5 x 10^ cm^"^) for model C30. Note that the box scale 
and level of grid are different for each panel. The spatial scale of each successive panel is 
different by a factor of four and thus the scale between panels (a) and (d) is different by a 
factor of 64. The magnetic field has an angle 6'^ ~ 30° in panel (a), where 6'^ is defined as 
the angle between the volume average magnetic field in the grid and the 2;-axis. Although the 
magnetic field fines in the thin disk run parallel to the disk surface {6b — 90°; Figure 13 &), 
the magnetic field fines outside the disk preserve the ambient direction of 9'^ ~ 30°. Closer 
to the cloud center, the magnetic field lines are twisted near the disk surface in the azimuthal 
direction, as seen in Figures 13 (c) and (d). As a result, the directions of the magnetic field 
are considerably different for different scales or densities even in the same cloud. 



3.3. Disk Formation Affected Both by Magnetic Field and Rotation 

In this subsection, we detail the evolution of group D, in which both the magnetic 
field and angular velocity are crucial for cloud evolution and disk formation. Thus, group 
D is located near the border between the magnetic-force- and rotation-dominant models, 
VLq/Bq ~ 0.39G'^/^c7^, and has the parameters a = 1 and u = 0.5. Group D has the same 
angular velocity as group C, but has a 10-times larger magnetic field. The clouds in this 
group have both a strong magnetic field and rapid rotation. 

The left panel of Figure 14 shows the directions of B, ft, and P for model D45. It 
can be seen that the direction of the magnetic field oscillates in the range 6b — 45° - 70°. 
Although the direction of the disk normal approaches the magnetic field direction, they are 
not completely aligned, as is the case for groups A and B. The direction of the rotation axis 
also oscillates considerably in the isothermal phase and is at = 63° at the core formation 
epoch. This value is almost the same as the angles of the magnetic field (^b = 51°) and the 
disk normal {9p — 49°), but the azimuthal coordinate 0n = 39° is very different from those 
of the magnetic field {(f) b = 283°) and the disk normal {(f) p = 275°). Thus, the direction 
of the rotation axis differs greatly from those of the magnetic field and disk normal. The 
direction of the disk normal is nearer to the magnetic field {if) bp = 6°) than the rotation 
axis {if)ap = 85°) at the core formation epoch, as shown in Table 2. Thus, the disk normal 
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seems to be parallel to the magnetic field for model D45 (Figure 14 right panel). On the 
other hand, the disk normal is nearer to the rotation axis {i/jcip = 6°) than the magnetic field 
{i^BP = 72°) in model D30. In model D60, the disk normal is close to both the magnetic 
field {ipBP = 29°) and the rotation axis {ipnp = 33°). In these clouds, the direction of the 
magnetic field, rotation axis, and disk normal oscillate in the isothermal phase. Whether a 
disk is aligned perpendicularly to the magnetic field or rotation motion depends on Oq. The 
cloud evolutions are influenced both by the magnetic fleld and the centrifugal force in group 
D. Therefore, we cannot clearly classify models in group D into either magnetic- force- or 
rotation-dominant models. 



4. Magnetic Flux— Spin Relation 

4.1. Amplification of the Magnetic Field and Angular Velocity 

The magnetic field strength and angular rotation speed increase as a cloud collapses. 
We have found in Paper II that the magnetic field strength normalized by the gas pressure 
and the angular velocity normalized by the free-fall timescale satisfy Equation (1) after 
a contracting disk forms in the isothermal phase for an aligned rotator model. In this 
subsection, we investigate the above relation for the case when the magnetic field is not 
necessarily parallel to the rotation axis (non-aligned rotator models). The evolution loci of 
the cores are plotted in Figure 15, where the horizontal and vertical axes were calculated 
using the central values of Pc, Be, and Qc- In Figure 15 a thick band indicating the equation 

r)2 r>2 

I = 1 (12) 

(0.2)2 ^^Qp^ ^ (0.36)2 8^c2p, ^ ' 

is also drawn, where the numerators of the left-hand side are defined as Be = {B^ ^ + By ^ + 
^2^)^/2 = (^s,c + ^y,c + ^z,c)^^^' where the suffix c indicates the values at the center. 

First, we consider the aligned rotator models (models with — 0) and compare them 
with those of Paper II. Comparing the solid lines (^o = 0° models) of Figure 15 with those 
of Figure 12 in Paper II, we can see the evolution locus is almost the same. That is, (1) the 
points move from the lower-left to the upper-right inside the B-Q relation line, whereas (2) 
those distributed outside the line move from the upper-right to the lower-left; (3) the slope 
of each evolution path is dlogQc/d^ogBc — 1 and (4) the evolution paths finally converge to 
Equation (12) for the isothermal phase. However, the evolution paths for models COO and 
DOO, which are located outside the B-Q relation line have slightly smaller angular velocities 
[Qc i'inGpe)-'^/^ ~ 0.1-0.15] than those of models C and D of Paper II [Qc (47rG'pc)"^/^ ^ 0.2]. 
This seems to be due to the fact that the initial cloud considered in this paper is more unstable 
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against gravity than that of Paper II; the ratio of thermal energy to gravitational energy 
is CKo = 0.168 in the present paper, while — 0.6-0.7 in Paper II. Although the clouds 
collapse vertically (vertical collapse) for the models C, D, COO, and DOO until the magnetic 
field and rotation satisfy the B-Q relation, the vertical collapse overshoots the B-Q relation 
line for initially unstable clouds in models COO and DOO. This difference can also be seen in 
comparing the density increase rate. The increase of the central density can be approximated 
by Pc ^ 5.1/[47rG'(t — t/)^] in model COO, where tf is the time at which the central density 
becomes infinite for the isothermal phase. This density increase rate is 1.75-times slower 
than that of the similarity solution (pc = 1.667/[47rG'(i — tf)^]; Larson 1969; Penston 1969). 
On the other hand, the density increase can be approximated by Pc ^ 6.2/[A7rG{t — tf^] for 
model C of Paper II, giving a density increase rate (6.2/5.1)^/^ ~ 1.1-times faster in model 
COO than in model C. This is a natural consequence of the lower ciq of model COO. This 
difference is also seen in the evolution of the oblateness. For example, the oblateness reaches 
^ob — 10 in model COO, while it reaches only Eoh — 4 in model C. Thus, a thinner disk is 
formed in model COO, which has a more unstable initial state. 

The evolution locus of model FOO moves towards the upper-right in the period Uc ^ 
lO^cm"^ in Figure 15. This indicates that the cloud collapses spherically in this period, 
as for models AOO and BOO, because the cloud is more unstable against the gravity than 
those of groups C and D. Group F has the same thermal energy as groups C and D, but 
has smaller rotational and magnetic energies, as shown in Table 1. Thus, in model FOO, 
the cloud collapses spherically in the early phase of the isothermal collapse. However, the 
collapse becomes anisotropic in the late phase of isothermal collapse because the magnetic 
force becomes effective. Although there are a few differences between the models, the con- 
vergence to the B-Q relation curve is evident, irrespective of the initial cloud shape and the 
distributions of the density, magnetic field, and angular velocity when the magnetic field is 
parallel to the rotation axis. This is natural because the B-Q relation is satisfied for the 
central part of the cloud and information on the outer part of the cloud is lost as the cloud 
collapses in the isothermal phase, as noted by Larson (1969). 

Next, we consider the non-aligned rotator models (^^o 7^ 0). In Figure 15, the dotted, 
broken, and dash-dotted lines show the evolution paths of models with 9 = 30°, 45°, and 60°, 
respectively. The points located inside the B-Q relation fine (non-aligned rotator models in 
groups A, B, and E) move from the lower-left to the upper-right, as for the aligned rotator 
models. These models have almost the same loci as the aligned rotator models and converge 
to the B-Q relation line. On the other hand, models located outside the B-Q relation fine 
(non-aligned rotator models in groups C, D, and F) show different evolutions than the aligned 
rotators. The non-aligned rotator models C30, C45, and C60 evolve towards the lower-right, 
then reverse their direction after they reach the B-Q relation line. Thus, the evolution of 
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the angular velocity in the non-aligned rotator models C30, C45, and C60 is the same as 
that of the aligned rotator model COO, while the evolution of the magnetic field is different. 
The magnetic field strength normalized by the gas pressure increases in non-aligned rotator 
models, but decreases in aligned rotator models. The evolution paths of non-aligned rotator 
models of group F (F40, F45, and F60) are toward the upper-left. Thus, the magnetic 
field strength normalized by the thermal pressure approaches the B-fl relation line in these 
models, while the angular velocity normalized by the free-fall timescale continues to increase 
in the isothermal collapse phase. The evolution paths for models D30, D45, and D60 oscillate 
in the Be {S7rclpc)~^^^-^c {^'^Gpc)~^^^ plane, moving away from the B-fl relation hue. 



4.2. Generalized Magnetic Flux— Spin Relation 

The growths of the magnetic field strength and angular velocity depend on the geometry 
of the collapse (vertical collapse to form a disk, spherical collapse, or lateral collapse in a 
disk). Figure 16 is similar to Figure 15, but with the magnetic field strength and the angular 
velocity parallel to the disk normal. That is, Bcp in the abscissa and in the ordinate are 
defined as 

Bcp^B- p, (13) 

p, (14) 

where B, Q,, and p represent the magnetic flux density vector, angular velocity vector, and 
the unit vector of the disk normal. The starting points of each locus are different even for 
models with the same a and u; but different (cf. AOO and A45), since the angle of the 
magnetic fleld at the initial stage and the formation epoch of the disk structure are dependent 
on ^0- Wc plotted the loci according to Equations (13) and (14) subsequent to the formation 
of a disk structure. The figure shows that all the evolution paths of models with the same 
a and uj (e.g. AOO and A45) move in the same direction irrespective of 6o, even though the 
starting points are different. 

Clouds for the models inside the B-fl relation line (support-deficient models: groups 
A, B, and E) evolve towards the upper-right, regardless of the initial angle 6o in Figure 16. 
Clouds having parameters inside the B-Q relation line collapse spherically until the magnetic 
field strength and angular velocity reach the B-fl relation line, as shown in §4.1. The 
clouds evolve isotropically (spherically), because the anisotropy caused by the magnetic or 
centrifugal force is not induced before the clouds reach the B-Q relation fine for small initial 
magnetic and rotational energies. Thus, the evolution direction of the non-aligned rotator 
models A45, B45, and E45 is the same as that of the non-ahgned rotator models AOO, BOO, 
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and EOO in Figure 15, though the angles between the magnetic field and rotation axis are 
different. This is because the anisotropy grows only shghtly during spherical collapse. 

For the support-sufficient models, the evolution paths of the non-ahgned rotator models 
C45, D45, and F45 are not the same as the aligned rotator models COO, DOO, and FOO in 
Figure 15. The cloud collapses vertically (vertical collapse) for groups C, D, and F, as shown 
in §3.2 and Paper II. In the case of the evolution of a weakly magnetized cloud rotating 
rapidly, in which the magnetic field is not parallel to the rotation axis and the magnetic 
field does not affect the cloud evolution, as seen in model C45, the cloud collapses along 
the rotation axis and the lateral collapse is suppressed by the centrifugal force. The cloud 
then forms a disk perpendicular to the rotation axis. The magnetic field and angular velocity 
parallel to the disk normal increase slightly for this collapse {flcp and B^p ~ constant). Thus, 
as the collapse proceeds, the evolution path is toward the lower-left, as shown in Figure 16 
and also seen for model COO. The magnetic field perpendicular to the disk normal (parallel 
to the disk) is amplified with the cloud collapse, when the cloud has a magnetic field that is 
not parallel to the disk normal at the initial stage. Including this component of the magnetic 
field, the evolution path is toward the lower-right, as shown in Figure 15 [the numerator of 
the second term of Equation (12) increases]. This is the reason why the normalized magnetic 
field strength in model C45 of Figure 15 increases in the isothermal collapse phase. The case 
for group F is similar to that of group C, however the cloud evolution is mainly controlled 
by the magnetic field for group F, not by the centrifugal force as in group C. Thus, the roles 
of the magnetic field and angular velocity are reversed. In group F clouds collapse along the 
magnetic field lines and disks are formed perpendicular to the magnetic field. The angular 
velocity perpendicular to the magnetic field or disk normal is then amplified and the angular 
velocity parallel to the magnetic field increases slightly. 

In Figure 16, the ordinate Qcp (47rG'pc)~^^^ indicates the ratio of the rotation and gravi- 
tational energies, while the abscissa Bcp (Sttc^Pc)"^''^ indicates the ratio of the magnetic and 
thermal energies. Since these ratios decrease in the aligned rotator models in proportion to 

— 1/2 

Pc for the vertical collapse phase, the thermal and gravitational energies catch up with 
the magnetic and rotational energies eventually. The cloud then reaches the B-Q relation 
line and the geometry of the collapse changes from vertical to lateral in the disk. Thus, a 
balance between the magnetic, rotational, thermal, and gravitational forces is achieved in a 
collapsing cloud. This type of evolution occurs in groups C, D, and F. 

When the evolution loci for groups A and B approach the B-Q relation line, the evolution 
depends on 9q even for the same a and cu. For groups A and B, both ratios (4'7rGpc)~^^^ 
and Be (Sttc^Pc)"^'^^ increase in proportion to pl^^ for a spherical collapse in which the cloud 
has small magnetic and rotational energies. Thus, the magnetic and rotational energies 
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become comparable to the gravitational and thermal energies during the contraction. The 
cloud then reaches the B-fl relation line and the geometry of the collapse changes from 
spherical to lateral. As a result, the anisotropy in 6 appears as the cloud reaches the B-Q 
relation line. Groups A and B form magnetic-force-dominant disks, in which the direction of 
the disk normal is controlled by the magnetic field. The magnetic field strength normalized 
by gas pressure does not increase or decrease after a cloud reaches the B-fl relation line in 
Figure 15. The angular velocity, however, can increase even after the cloud has reached the 
B-^} relation line in Figure 15, because the rotation axis is not parallel to the disk normal. 
Differences in the position of the end point in models AOO, A30, A45, and A60 are caused 
by this mechanism. Magnetic braking is also effective in these models. For these reasons, 
the evolution paths begin to diverge as they approach the B-D, relation line. 

We have plotted the normalized magnetic field strength and angular velocity at the 
initial stage for models WF, MF, and SF of MT04 as diamonds in Figure 16. MT04 shows 
that disks are formed perpendicular to the local magnetic field in all non-aligned rotator 
models of MT04. This is natural because the models WF, MF, and SF are distributed in 
the magnetic-force-dominant region in Figure 16. 

In summary, the geometry of the collapse determines the amplification of the mag- 
netic field and angular velocity. The gas cloud in the support- deficient region amplifies the 
magnetic field strength and angular velocity during the contraction, regardless of the initial 
angle 6q. In these models, aligned and non-aligned rotators evolve similarly. On the other 
hand, for the support-sufficient models, aligned and non-aligned rotators evolve differently. 
In the rotation-dominated models, the magnetic field perpendicular to the rotation axis is 
amplified. This plays a role as a non-axisymmetric perturbation in forming a bar or spiral 
structure. Even in non-aligned rotator models, the generalized magnetic flux-spin relation 
holds in contracting disks formed in the isothermal regime. 

4.3. Disk Formation by Magnetic Field or Rotation 

We have shown that a cloud forms either a magnetic-force-dominant or a rotation- 
dominant disk according to its initial conditions and that cloud evolution can be well under- 
stood and classifled using the generalized B-D, relation. In this subsection, we show how we 
can specify the parameter regions where a disk is formed under the influence of either the 
magnetic fleld or under the influence of rotation. The magnetic field, rotation axis, and disk 
normal have different directions when the magnetic field is not parallel to the rotation axis 
at the initial stage, although they have are identical for an aligned rotator. We can assess the 
dominant force for disk formation using the evolution loci of the direction of the magnetic 
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field, rotation axis, and disk normal. When a disk is formed perpendicular to the magnetic 
field, the disk normal moves in association with a small-scale magnetic field in magnetic- 
force-dominant models (Figure 3). A similar evolution is seen in rotation-dominated models. 
We compare the evolution of 16 models with different a and io but the same Q^. We choose 
the initial angle between the magnetic field and rotation axis as = 45°, because the cloud 
evolution does not depend on ^o, as shown in §3.1 and §3.2. The angles (6*5, 6*^, 6'p), (0_b, 
0p)) and (t/'so, ^/'sp, ipfip)', the dominant force for forming a disk {B or Q); and the 
axis-ratio (ear) at the end of the isothermal phase are all listed in Table 2. The dominant 
force for disk formation is determined by the loci of the magnetic field, rotation axis, and 
disk normal for the isothermal phase. For example, the locus of the disk normal moves 
together with that of the magnetic field in magnetic-force-dominant models. 

The shapes of the clouds at the core formation epoch are shown in Figure 17. In this 
figure, each panel is positioned based on the initial magnetic field strength and angular 
velocity. The figure shows that the disk normals are almost parallel to the 2;-axis in the 
upper-left region, while they are parallel to the magnetic field in the lower-right region. This 
is natural since the cloud is initially rotating rapidly and magnetized weakly in the upper-left 
region, while it rotates slowly and is magnetized strongly in the lower-right region. 

In order to compare the cloud evolutions with different initial angular velocities we focus 
on four models with a = 0.01 and different uj = 0.3, 0.1, 0.03, and 0.01 [models (b), (f), (j), 
and (n)], shown in Figure 17. These models are aligned in the second column of Figure 17. 
The disk normal is oriented with the ^-direction, and the magnetic field lines arc inclined 
from the z-a.xis in the models with large uj [(b) and (f)]. Each model has a small angle 
between the rotation axis and the z-axis [9q = 2° (b), 3° (f), 2° (j), and 5° (n) ]. Thus, it is 
shown that the cloud evolves maintaining the direction of the initial rotation axis. On the 
other hand, the angle between the magnetic field and the 2;-axis increases with increasing 
initial angular velocity u [9b — 86° (b), 83° (f), 54° (j), and 46° (n)]. The magnetic field in 
models (b) and (f) is almost perpendicular to the rotation axis and the disk normal {i^BQ, 
ifjBp — 90°). It appears that the disk is formed by the effect of rotation in models (b) and 
(f), because the angle between the rotation axis and disk normal is small [V'f^p = 2° (b) and 
6° (f)]. On the other hand, the disk seems to be formed by the magnetic force in models (j) 
and (n), because the angles iI^bp [i^BP — 25° (j) and 4° (n)] are smaller than those of ■0np 
[i^vip — 28° (j) and 39° (n)]. The axis-ratio increases with increasing u [sar = 0.51 (b), 0.18 
(f), 0.91 X 10~^ (j), and 6 x 10~^ (n)], because the disk forms earlier in a model with larger 

UJ. 

Next, we focus on four models with the same u — 0.1 but different a — 0.001, 0.01, 0.1, 
and 1 [models (e), (f), (g), and (h)], in order to compare cloud evolution with different initial 
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magnetic field strengths. These models are ahgned in the second row of Figure 17. The disk 
normals are considerably inclined from the z-axis in the models with strong magnetic field, 
(g) and (h). This inclination indicates that the disk is formed by the effect of the magnetic 
force. The disk is perpendicular to the rotation axis in model (e) [{ipBP, i^np) = (83°, 1°)], 
while the disk is perpendicular to the local magnetic field rather than the rotation axis in 
model (h) [{^BP, M = (2°, 11°)]. 

The angle between the rotation axis and the disk normal {ipnp) is smaller in model (h) 
(■^np = 11°) than in model (g) {ipnp = 48°). This seems to be due to the fact that the 
angular momentum perpendicular to the magnetic field in models with a > 0.1 is effectively 
removed by the magnetic braking process, and thus the direction of rotation tends to incline 
to the magnetic field and the disk normal. As a result, the rotation axis is considerably 
inclined from the initial direction in models with a strong magnetic field [9q = (e), 3 
(f), 16 (g), and 32 (h)]. The magnetic field is parallel to the disk in models with small a 
[e.g. i/jbp = 83 (e)], while the magnetic field maintains its initial direction in models with 
large a [e.g. 9b — 43 (h)]. The azimuthal directions of the magnetic field (0b) and disk 
normal {4>p) coincide in models (e), (f), (g), and (h). In these models, the gas contracts 
along the magnetic field line onto the disk mid-plane, then a non-axisymmetric structure 
(i.e. bar structure) is formed perpendicular to the magnetic field, as discussed in §3.2. The 
non-axisymmetry tends to increase with the initial magnetic field strength [see models (e), 
(f), and (g)], except for model (h). 

The role of the magnetic field in magnetic-forcc-dominant models is similar to that of the 
centrifugal force in rotation-dominant models. The disk orientation is essentially determined 
by the direction of the dominant force. However, there is at least one quantitatively different 
point between the two types of models. Namely, the angular momentum is transferred by the 
magnetic braking in the magnetic-dominant models. The dominant force (B or Q) for disk 
formation is summarized in Table 2 and Figure 17. The shadowed region in the lower-right 
part of Figure 17 indicates disks formed by the Lorcntz force, while the upper-left region 
indicates disks formed by the centrifugal force. A broken line between these two indicates 
the border between the magnetic-force-dominant and rotation-dominant disks, which is well 
fitted by 

^^0.39G'/\-\ (15) 

-Do 

similar to the result for aligned rotators. 

Cloud evolution can be classified into four patterns using the generalized B-fl relation 
curve 

' +77V^J^^ = 1. (16) 



(0.2)2 47rG'pc (0.36)2 87rc2pe 
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and Equation (15): (i) support-deficient, rotation-dominant models [inside the B-D, relation 
line and above Equation (15)], (ii) support-deficient, magnetic-force-dominant models [inside 
the B-Q relation line and below Equation (15)], (iii) support-sufficient, rotation-dominant 
models [outside the B-Q relation line and above Equation (15)], and (iv) support-sufficient, 
magnetic-force-dominant models [outside the B-Q relation line and below Equation (15)]. In 
the models of class (i) [models (e), (f), (g), and (i)], the cloud collapses slowly, maintaining 
spherical symmetry, and then a disk forms due to the rotation. In this type of evolution, 
the magnetic field hardly changes its initial direction. On the other hand, in the models 
of class (iii) [models (a), (b), and (c)], the cloud collapses along the rotation axis owing to 
the strong centrifugal force, and then a thin disk is formed in the early isothermal collapse 
phase. The magnetic field lines run along the disk plane, because the magnetic field lines 
are compressed together with a cloud in these models. In the rotation-dominant models of 
classes (i) and (iii), the rotation axis maintains its initial direction because the magnetic 
braking is not so effective. On the other hand, in the magnetic-forcc-dominant models of 
classes (ii) and (iv), the rotation axis is inclined from the initial direction because the clouds 
with strong magnetic fields experience effective magnetic braking. The inclination of the 
rotation axis, 9q in class (iv) is greater than that in class (ii). The direction of the magnetic 
field, however, tends to maintain its initial direction relative to the z-axis in classes (ii) and 
(iv) for a strong magnetic tension force. 



5. Discussion 

5.1. Fragmentation of a Magnetized Rotating Cloud 

Fragmentation is considered to be one of the mechanisms producing binary and multiple 
stars. We investigated fragmentation of a rotating magnetized cloud in Paper III for the case 
of Bq W Oq. We found fragmentation only in rotation-dominant clouds. This indicates that 
fragmentation is suppressed by the magnetic field, similar to the findings of Hosking & 
Whitworth (2004) and Ziegler (2005). Fragmentation occurs via a global bar or ring mode 
and depends on the initial amplitude of the non-axisymmetric perturbation in support- 
sufficient rotation-dominant clouds. That is, when the non-axisymmetric structure barely 
grows in the isothermal phase and the rotation rate reaches uj ^ 0.2 at the core formation 
epoch, an adiabatic core fragments via a ring (ring fragmentation). On the other hand, when 
the core is deformed to an elongated bar at the core formation epoch, the bar fragments 
into several cores (bar fragmentation). The parameter study in Paper III shows that ring 
fragmentation is seen in rot at ion- dominant models in the adiabatic phase [classes (i) and 
(iii)] and bar fragmentation is observed only in support-sufficient rotation-dominant models 
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[class (iii)]. Support-deficient and support-sufficient magnetic-dominant models [classes (ii) 
and (iv)] evolve into a single dense core without fragmentation owing to effective magnetic 
braking and a slow rotation. 

The results obtained in Paper III show that fragmentation patterns are dependent on 
the growth of a non-axisymmetric structure. The non-axisymmetric perturbation begins to 
grow after a thin disk is formed. The amplitude of the non-axisymmetric mode barely grows 
in support-deficient rotation-dominant models [class (i)], because the disk forms slowly, as 
shown in §3.1. On the other hand, the non-axisymmetric structure grows sufficiently in 
support-sufficient rotation-dominant models [class (iii)], because the disk forms promptly. 
In class (iii), the patterns of fragmentation are dependent on the initial amplitude of the 
non-axisymmetric perturbation. Thus, a cloud fragments through a ring when the cloud 
has a small amount of initial non-axisymmetric perturbations, whereas the cloud fragments 
through a bar when it has a sufficient amount of initial non-axisymmetric perturbations. 
These results apply to aligned rotator clouds. Below we extend the study to the evolution 
of non-aligned rotator models. 

We did not add any explicit non-axisymmetric perturbations to the initial state of non- 
aligned rotator models. We did not find any rings for the rotation-dominant models of 
non-aligned rotators. In the rotation-dominant models, non-axisymmetry arises from the 
magnetic force in the case of Bq j]' Qq. Since this non-axisymmetric perturbation from the 
magnetic force grows sufficiently in the isothermal phase, bar fragmentation must occur in 
non-aligned rotator models. 

As detailed previously, the axis-ratio (^ar), listed in Table 2, arises from an anisotropic 
force due to the magnetic field or rotation. Model (c) has the greatest axis-ratio ~ 2.2 
of all models. The initial state of this cloud is outside the B-Q relation line, and has the 
strongest magnetic field in rotation-dominant models (Figure 17). In this cloud, a disk is 
formed perpendicular to the rotation axis and the magnetic field parallel to the disk surface 
is greatly amplified. This induces a bar structure along the magnetic field line on the disk 
mid-plane, as shown in §3.2 (e.g. model C45). Although this bar does not fragment in this 
study, such a bar structure suggests the possibility of fragmentation in the adiabatic phase 
(Paper III) if the calculation is continued. Even if the bar does not fragment in the case of 
there being no explicit non-axisymmetric perturbations, we expect bar fragmentation if an 
initial explicit perturbation is added. We confirmed that bar fragmentation occurs in model 
(c) when we added a 10% non-axisymmetric density perturbation to the initial state. As 
a result, the anisotropy arising from magnetic and centrifugal forces in non-aligned rotator 
models promotes bar fragmentation and suppresses ring fragmentation. 
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5.2. Compcirison with Previous Works 

Several studies have examined the gravitational collapse of molecular cloud cores and 
star formation using three-dimensional MHD calculations (Dorfi 1982, 1989; Boss 2002; Hosk- 
ing & Whitworth 2004; Ziegler 2005; Banerjee & Pudritz 2006, MT04; Papers I, II and III). 
Except for those of Dorfi (1982, 1989) and MT04, these studies assume that the magnetic 
field lines are initially parallel to the rotation axis. The evolutions of non-ahgned rotator 
models are almost the same as those of aligned rotator models for support-deficient models, 
while the evolutions are completely different for support-sufficient models. For example, the 
direction of the magnetic field continues to move away from the rotation axis in a rapidly 
rotating cloud. The directions of the magnetic field and rotation are completely different af- 
ter disk formation in this case. Comparing magnetic fields parallel and perpendicular to the 
rotation vector shows that a perpendicular magnetic field transfers angular momentum more 
effectively than a parallel field (Mouschovias & Paleologou 1979). In other words, magnetic 
braking is more effective in non-aligned rotator models than in aligned rotator models. This 
strong magnetic braking may be the solution of the angular momentum problem, in which 
the specific angular momentum of a parent cloud is much larger than that of a new-born 
star. 



5.3. Comparison with Observation 

The magnetic field strengths and directions have been observed for many clouds. It 
is believed that there is no correlation between the direction of the magnetic field and the 
large-scale cloud shape (Goodman et al. 1993; Tamura et al. 1995; Word- Thompson 2000). 
Recently, the directions of the magnetic field have been observed in both large (cloud) and 
small (prestellar core) scales in the same target. The small-scale polarization pattern of the 
W51 molecular cloud observed by BIMA (Berkeley-Illinois-Maryland Association; Lai et al. 
2001) coincides with the large-scale polarization pattern observed by SCUBA (Chrysostomou 
et al. 2002). Also, the directions of the magnetic field were found to be the same in both large- 
and small-scale observations of the DR21 cloud (Lai et al. 2003). However, some observations 
have given contrasting findings. Although the average polarization angle in the MMS6 core in 
the OMC-3 region of the Orion A cloud (Matthews et al. 2005) coincides with the large-scale 
polarization angle observed by SCUBA (Houde et al. 2004), the polarization angle changes 
systematically across the core. An observation of the Barnard 1 cloud in Perseus reveals that 
three of the four cores exhibit different mean field directions than that of the ambient cloud 
(Matthews & Wilson 2002; Matthews et al. 2005). These trends of OMC-3 and the Barnard 
1 cloud agree well with the results for non-aligned rotator models in group C in §3.1. The 
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direction of a small-scale magnetic field can be different from the large-scale field in models 
C30, C45, and C60 (support-sufficient, rotation-dominant models). As shown in Figure 13 
for model C30, although the magnetic field maintains its initial direction Ob — 30° outside the 
high-density core, inside the high-density region the magnetic field lines are perpendicular 
to the rotation {z-) axis {9b — 90°). Thus, the direction of the magnetic field varies for 
different spatial scales. In this model, the angle between the large-scale and small-scale 
magnetic fields is 85°. On the other hand, the magnetic field fines have the same direction in 
both large- and small-scales in W51 and DR21. These clouds correspond to groups A, B, and 
E, in which the magnetic field hardly changes its direction in the isothermal phase. These 
clouds are expected to have a slow rotation rate. These findings show that the direction of 
the magnetic field can change only in support-sufficient rotation-dominant clouds. 

Recently, an hour-glass-shaped magnetic field has been found in a dynamically con- 
tracting core around the binary protostellar system NGC 1333 IRAS 4A (Girart et al. 2006). 
Two outfiows were also observed in this region and associated with each protostar of the 
protobinary system (Choi 2005). However, the direction of the magnetic field docs not coin- 
cide with the outflow axis (Girart et al. 2006). From our previous study, fragmentation (or 
binary formation) appears only in rotation-dominated clouds (Paper III). In these clouds, 
the magnetic fleld tends to be aligned along the direction perpendicular to the rotation axis, 
as shown in §3.2.1, and therefore the direction of the magnetic field in a dense core does not 
coincide with that of the large-scale field. Thus, the observed miss- aligned outfiow indicates 
that a binary is being formed from a rotation-dominated cloud, because the outfiows are 
driven along the local magnetic field (MT04). 

Our numerical calculations were carried out with a Fujitsu VPP5000 at the Astronomical 
Data Analysis Center of the National Astronomical Observatory of Japan. This work was 
supported partially by Grants- in- Aid from MEXT (16077202 [MM], 17340059 [TM, KT], 
15340062, 14540233 [KT], 16740115 [TM]). 
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Table 1: Parameters and initial conditions for typical models 



Group Model a u 6^ f ao 7o "o" ^o'' ^o^ M'^ 

A AOO, A30, A45, A60 0.01 0.01 (0°, 30°, 45°, 60°) 1.05 0.70 3.29x10-^ 1.34x10-2 5.25 3.23 1.38 6.41 
B BOO, B30,B45,B60 0.1 0.01 (0°, 30°, 45°, 60°) 1.05 0.70 3.29x10-^ 0.134 5.2510.21.38 6.41 
C COO, C30, C45, C60 0.0 1 0.5 (0°, 30°, 45°, 60°) 5.0 0.168 0.823 3.22x10"^ 25 6.59 141 26.7 
D DOO, D30, D45, D60 1 0.5 (0°, 30°, 45°, 60°) 5.0 0.168 0.823 0.32 25 65.9 141 26.7 

E E00,E30,E45,E60 10-^0.05(0°, 30°, 45°, 60°) 5.0 0.168 0.82x10-^^3.22x10-^ 25 2.0814.126.7 
F FOO, F30, F45, F60 1 0.05 (0°, 30°, 45°, 60°) 5.0 0.168 0.82x10"^ 0.32 25 65.914.126.7 
SF^ 3.04 0.14 (00°, 45°, 90°) 1.68 0.5 002 2^88 2.6137.17.116.13 

MF^ 0.76 0.14 (45°, 70°, 80°) 1.68 0.5 0.02 0.72 2.61 18.6 7.116.13 

WF^ 0.12 0.14 (00°, 45°, 90°) 1.68 0.5 0.02 0.12 2.617.42 7.116.13 

« no(10^ X cm.-^), ^ Bq (/nG), ^ (10"'^ yr^^), M (M©) , " models calculated by Matsumoto & 

Tomisaka (2004) 



Table 2: Calculation results at the core formation epoch 



Model 


a 




^^0 


/ 


(Ob; dn, Op) 


(0s, 0p) 


(^Bo, V'sp, -ipnp) 








AOO 


0.01 


0.01 


00 


1.05 


(0, 0, 0) 


(0, 0, 0) 


(0, 0, 0) 


— 







A30 


0.01 


0.01 


30 


1.05 


(31, 1, 30) 


(24, 22, 20) 


(30, 2, 29) 


B 


1.2 X 10- 


-3 


A45 


0.01 


0.01 


45 


1.05 


(46, 2, 44) 


(24, 23, 20) 


(44, 4, 42) 


B 


7.1 X 10" 


-3 


A60 


0.01 


0.01 


60 


1.05 


(61, 2, 58) 


(24, 11, 20) 


(59, 4, 56) 


B 


1.5 X 10" 


-2 


BOO 


0.1 


0.01 


00 


1.05 


(0, 0, 0) 


(0, 0, 0) 


(0, 0, 0) 


— 







B30 


0.1 


0.01 


30 


1.05 


(33, 32, 33) 


(22, 66, 14) 


(24, 4, 27) 


B 


8.0 X 10- 


-3 


B45 


0.1 


0.01 


45 


1.05 


(46. 7. 45) 


(23, 14, 15) 


(39, 5, 38) 


B 


2.5 X 10- 


-3 


B60 


0.1 


0.01 


60 


1.05 


(58, 52, 57) 


(24, 298, 15) 


(68, 8, 61) 


B 


1.3 X 10- 


-2 


COO 


0.01 


0.5 


00 


5 


(0, 0, 0) 


(0, 0, 0) 


(0, 0, 0) 


— 







C30 


0.01 


0.5 


30 


5 


(88, 1, 2) 


(302, 273, 96) 


(88, 90, 3) 


n 


0.43 




C45 


0.01 


0.5 


45 


5 


(90. 1. 2) 


(302, 263, 86) 


(88, 89, 3) 


n 


0.61 




C60 


0.01 


0.5 


60 


5 


(89, 1, 2) 


(302, 257, 84) 


(89, 89, 3) 


n 


0.68 




DOO 


1 


0.5 


00 


5 


(0, 0, 0) 


(0, 0, 0) 


(0, 0, 0) 


— 







D30 


1 


0.5 


30 


5 


(60, 13, 14) 


(243, 8, 34) 


(67, 72, 6) 


n 


0.19 




D45 


1 


0.5 


45 


5 


(51, 63, 49) 


(283, 39, 275) 


(88, 6, 85) 


B 


0.21 




D60 


1 


0.5 


60 


5 


(35, 26, 51) 


(337, 51, 13) 


(36, 29, 33) 


B 


0.32 




EOO 


0.001 


0.05 


00 


5 


(0, 0, 0) 


(0, 0, 0) 


(0, 0, 0) 


— 







E30 


0.001 


0.05 


30 


5 


(43, 1, 1) 


(50, 80. 38) 


(43, 42, 1) 




3.4 X 10- 


-3 


E45 


0.001 


0.05 


45 


5 


(59, 1, 1) 


(50, 110, 47) 


(58, 58, 1) 




6.8 X 10- 


-3 


E60 


0.001 


0.05 


60 


5 


(71, 1, 1) 


(90, 142, 49) 


(71, 70, 1) 




1.3 X 10- 


-2 


FOO 


1 


0.05 


00 


5 


(0, 0, 0) 


(0. 0. 0) 


(0, 0, 0) 


— 







F30 


1 


0.05 


30 


5 


(35, 9, 30) 


(51, 13, 24) 


(28, 15, 21) 


B 


0.21 




F45 


1 


0.05 


45 


5 


(51. 10. 44) 


(51. 354. 24) 


(46. 20. 35) 


B 


0.40 




F60 


1 


0.05 


60 


5 


(65. 12. 58) 


(5L 336. 21) 


(62, 21, 50) 


B 


0.10 




(a) 


0.001 


0.3 


45 


1.68 


(88, 0, 0) 


(189, 308, 224) 


(88, 88, 0) 




5.9 X 10- 


-2 


(b) 


0.01 


0.3 


45 


1.68 


(86, 2, 0) 


(294, 269, 87) 


(85, 87, 2) 




0.51 




(c) 


0.1 


0.3 


45 


1.68 


(82, 15, 8) 


(226, 356, 167) 


(89, 78, 23) 


Q 


2.2 




(d) 


1 


0.3 


45 


1.68 


(48, 46, 42) 


(256, 32, 231) 


(85, 18, 87) 


B 


0.17 




(e) 


0.001 


0.1 


45 


1.68 


(83, 0, 0) 


(120, 42, 127) 


(83, 83, 1) 


Q 


1.7 X 10- 


-2 


(f) 


0.01 


0.1 


45 


1.68 


(83, 3, 4) 


(125, 4, 126) 


(85, 79, 6) 


n 


0.18 




(g) 


0.1 


0.1 


45 


1.68 


(78, 16, 37) 


(146, 347, 109) 


(87, 50, 48) 


n 


0.44 




(h) 


1 


0.1 


45 


1.68 


(43, 32, 42) 


(63, 68, 60) 


(12, 2, 11) 


B 


5.2 X 10- 


-2 


(i) 


0.001 


0.03 


45 


1.68 


(53, 0. 4) 


(51, 81, 49) 


(53, 49, 4) 




2.1 X 10- 


ii 

~ A 


(j) 


0.01 


0.03 


45 


1.68 


(54, 2, 30) 


(53, 25, 45) 


(52, 25, 28) 


B 


9.1 X 10- 


-2 


(k) 


0.1 


0.03 


45 


1.68 


(58, 5, 45) 


(70, 342, 35) 


(58, 30, 40) 


B 


6.0 X 10- 


-2 


(1) 


1 


0.03 


45 


1.68 


(45. 26. 45) 


(20, 16, 19) 


(19, 1, 19) 


B 


5.3 X 10- 


-2 


(m) 


0.001 


0.01 


45 


1.68 


(46, 0, 29) 


(18. 2. 16) 


(46, 17, 29) 


B 


1.1 X 10- 


-2 


(n) 


0.01 


0.01 


45 


1.68 


(46, 5, 43) 


(19, 21, 16) 


(41, 4, 39) 


B 


6.0 X 10- 


-3 


(o) 


0.1 


0.01 


45 


1.68 


(45, 8, 45) 


(29, 4, 12) 


(38, 12, 37) 


B 


5.8 X 10" 


-3 


(P) 


1 


0.01 


45 


1.68 


(45, 29, 45) 


(6, 5, 6) 


(16, 0, 16) 


B 


2.8 X 10- 


-2 



"The dominant forces for disk formation. They are determined by the loci of the magnetic field, the rotation 
axis, and disk normal for the isothermal phase. 



- 31 - 




Fig. 1. — Time sequence of model A45. The structure of the high-density region (n > O.lric ; 
isosurface), density contours (contour hues), velocity vectors (arrows), and magnetic field 
lines (streamlines) are plotted in panels (a)-(d). The green colored disk indicates the region 
of n > 1/100 nc on the mid-plane parallel to the disk-like structure (perpendicular to the 
disk normal). The panels show the stages (a) ra^ = 5 x lO^cm"^ (initial stage), (b) = 
5.9 X lO^cm"^, (c) ric = 5.5 x 10^°cm~^, and (d) = 2.5 x lO^'^cm"^. The level of the finest 
grid (/), elapsed time from the beginning (t), and central number density {ric) are hsted at 
the top of each panel. The size of the grid is also shown. 
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Fig. 2. — Oblateness (upper panel) and axis-ratio (lower panel) against central density for 
group A (models AOO, A30, A45, and A60). 
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Fig. 3. — Loci of directions of the magnetic field B (dotted line), rotation axis fl (broken 
line), and disk normal P (solid line) for model A45. The symbols o, o. A, and * denote the 
stages ric = 10® cm~^, 5 x 10^'^ cm~^, 10^^ cm~^, and the final epoch, respectively. 
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Fig. 4. — Angles 9b (broken line), 6q_ (solid line), and 9p (dotted line) for models AOO, A30, 
A45, and A60, plotted against the central number density (ric). The angles 0b (thin broken 
line), (thin solid line), and 0p (thin dotted line) for model A45 are also plotted. The left 
and right ordinates indicate the angles 9 and 0, respectively. The abscissa represents the 
central number density. 
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320 AU 



Fig. 5. — (Left) Loci of magnetic field, rotation axis, and disk normal for model B45. The 
symbols o indicate the angles at the core formation epoch. (Right) Structure of high density 
region (n > OAric, isodensity surface), density contours (contour lines), velocity vectors 
(arrows), and magnetic field line (streamline) for Uc = 5 x 10^° cm~^. 
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Fig. 6. — As Figure 5 for model F45 at Uc = 1.3 x 10^^ cm ^ (right). 




Fig. 7. — As Figure 5 for model E45 at Uc = 1.8 x 10^^ cm ^ (right). 
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Fig. 8. — Snapshots of model C45. The structure of the high-density region [n > O.lric] 
isodensity surface), density contours (contour hues), velocity vectors (arrows), and magnetic 
field lines (streamlines) are plotted in panels (a)-(d). The upper panels show the view along 
the y-axis (edge-on view). The lower panels show the view along the z-axis (face-on view). 
The green colored disk in panels (c) and (d) indicates the region of n > (1/100) ric on the mid- 
plane parallel to the disk-like structure. The panels show the stages (a) ric = 6.3 x lO^cm"^, 
(b) ric = 6 X lO^cm^^, (c) Uc = 2.3 x lO^^cm"^, and (d) Uc = 6.9 x lO^^'cm"^. The level of 
the finest grid (/), elapsed time from the beginning (t), and central number density {uc) are 
listed at the top of each panel. The size of the grid is also shown. 
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Fig. 10. — As Figure 5 (left) for model C45. The lower left inset is an enlarged view of the 
center. 
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Fig. 11. — As Figure 1 for COO at the core formation epoch. 




Fig. 12.— As Figure 4 for group C (models COO, C30, C45, and C60). 
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Fig. 13. — Views of model C30. Each panel shows the same epoch, t = 1.52 x 10^ yr 
{ric = 1.5 X lO^cm-^), but for different box scales, (a) 9.2 x 10^ AU, (b) 2.3 x 10^ AU, 
(c) 5.8 X 10^ AU, and (d) 1.4 x 10'^ AU. Isodensity surfaces are drawn at the densities (a) 
ric > 1.3 X lO^cm-3, (b) > 5 x 10^ (c) > 5 x 10^ and (d) > 5 x 10^ cm'^. The 
density distribution on the y = plane is projected on to the boundary by a false color plot. 
The other representatives are the same as those in Figure 1. 
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Fig. 14. — As Figure 5 for model D45 at rzc = 1.8 x 10^^ cm ^ (right). 
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Fig. 15. — Evolution of magnetic flux density and angular velocity at the cloud center. The 
X-axis indicates the square root of the magnetic pressure [i?c /(Svr)^/^] divided by the square 
root of the thermal pressure [(c^pc)^^^]- The y-axis represents the angular speed (f^c) divided 
by the free-fall rate [(47rGpc)^'^^]- Here, we use the definitions Be = {Bl^ + B^^ + 5^ jV2 



and f2r 



[a 



2 

x,c 



1^/^, where the suffix c indicates the value at the center. The 



upper X-axis indicates the value of a. The symbols *, A, o, □, and + represent the magnetic 
field and the angular velocity at ric = 5 x lO^cm"^ (initial state), lO^cm"^, lO^cm"^, 
lO^cm"^, and 5 x 10^°cm~^, respectively. Each line denotes the evolution path from the 
initial state (nc,o = 5 x lO^cm"^) to the end of the isothermal phase {ric = 5 x 10^°cm~'^). 
The characters A, B, C, D, and E denote the group names shown in Table 1. The different 
lines indicate different initial angles, 9q = 0° (solid line), 30° (dotted line), 45° (broken line), 
and 60° (dash dotted line). The thick grey band denotes the magnetic flux-spin relation 
nl/[{0.2f X ATiGpc] + 52/[(0.36)2 X Snclpd = 1 [see Equation (12)]. 
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Fig. 16. — Evolution of magnetic flux density and angular velocity parallel to the disk normal. 
The axes and symbols are the same as for Figure 15. Here, B^p and Qcp are the magnetic 
field and angular velocity parallel to the disk normal [Equations (13) and (14)]. Models with 
= 0° and 45° are plotted. The symbols o show the initial states of the models WF, MF, 
and SF calculated in MT04 (see Table 1). 
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Fig. 17. — Snapshots of a disk and magnetic field lines at ric = 5 x 10^*^ cm~^. The snapshots 
are displayed according to the initial values of flc and Be- The axes and curve are the 
same as that of Figure 15. The structure of the high-density region (n > O.lric; isodensity 
surface), density contours (contour lines), velocity vectors (arrows), and magnetic field hues 
(streamlines) are plotted in the panels (a)-(p). The shaded part indicates the region in 
which a magnetic-dominant disk is formed. The broken line denotes the border between the 
rotation- and the magnetic-force-dominant disks, Vt/B = 0.396"^/^ c^^ [Equation (15)]. 



